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We analyze measures of the azimuthal asymmetry, in particular the participant harmonic mo- 
ments, £*, in a variety of Glauber-like models for the early stage of collisions at RHIC. Quan- 
titative comparisons indicate substantial model dependence for £*, reflecting different effective 
number of sources, while the dependence of the scaled standard deviation (7(e*)/e* on the par- 
ticular Glauber model is weak. For all the considered models the values of (7(e* )/ e* range from 
^ 0.5 for the central collisions to ~ 0.3-0.4 for peripheral collisions. These values, dominated by 
statistics, change only by 10-15% from model to model. For central collisions and in the absence 
of correlations between the location of sources we obtain through the use of the central limit 



theorem the simple analytic formula a{e*)/ e* {b = 0) ~ y/A/n— 1 ~ 0.52, independent on the 
collision energy, mass number, or the number of sources. We investigate the shape-fluctuation 
effects for jet quenching and find they are important only for very central events. Finally, we list 
some remarks and predictions from smooth hydrodynamics on higher flow coefficients and their 
fluctuations, in particular (7(v4)/v4 ~2a{v2)/ vi- 
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Figure 1: A typical gold-gold collision in the x — y plane at b = 6 fm. Left; wounded nucleons. Red 
and black circles indicate nucleons from the two colliding nuclei. Right: the centers of mass of pairs of 
nucleons undergoing binary collisions. The straight lines indicate the twisted and shifted principal axes of 
the quadrupole moment, while the blue dots show the center of mass of the system. 



1. Introduction 

This talk is based on Ref. where more technical details may be found. 

Event-by-event hydrodynamic studies ||] of relativistic heavy-ion collisions have revealed 
that fluctuations of the initial shape of the system formed in the early stage of the reaction lead 
to quantitatively relevant effects for signatures of the azimuthal asymmetry [Q, |5|, ^, ^ ^]. These 
effects are also important for experimental measurements of the elliptic flow |l0[ [11], 12]. In 
this talk we report our investigation of this phenomenon in the framework of various Glauber-like 
approaches describing the deposition of energy in the system in the early stages of the collision. 
Our study focuses on both understanding of the statistical nature of the results, as well as on com- 
parisons of various models. 

Figure |] illustrates the two effects of the shape fluctuations due to the finite number of sources: 
the shift of the center-of-mass and the rotation of the the quadrupole principal axes. Statistical 
analyses may be performed in the reference frame fixed by the reaction plane (we call itfixed-axes, 
a.k.a. standard), or for each event in the frame defined by the twisted and shifted principal axes 
(we call it the variable-axes frame, a.k.a. participant). In the fixed-axes frame the two-dimensional 
probability distribution of sources can be Fourier-expanded as 



/(P, </>) = /o(p) + 2/2(p) cos(2(^>) + 2/4(p) cos(4(^)) + . . . ^ 



(1.1) 



where the transverse radius p = \Jx^ is measured from the center of the geometric intersection 
of the two nuclei. One also introduces 

\l%pfi{p)p^dp 



On the other hand, in the variable-axes frame we have the distribution 

r (p, 0) = /o(p) + 2/2*(p) cos(2(/) - 20*) + 2/4* (p) cos(40 - 4r ) + . . . , 
where <p* denotes the rotation angle of the principal axes in each event. Correspondingly, 

JlKpflipJp^ 
' j27ipfo{p)p^dp 



(1.2) 



(1.3) 



(1.4) 
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(/o = /q ). The quadrupole parameters are denoted without the subscript as £ = £2 and £* = £|. 
2. The toy problem 

Consider the one-dimensional problem where uncorrelated particles are randomly generated 
from a distribution in the azimuthal angle <p containing the monopole and quadrupole moments, 

/(0) = l+2£cos(20), eeh^^l- ^^-^^ 
Obviously, the distribution has only two non-zero fixed-axes moments, 

/o = 2^ j|'''#/(0) = l, f2 = ^fJd^cos{2<^)m=e. (1.1) 



We generate n particles according to the distribution (2.1) in each event, and subsequently carry 



out the averaging over the events, denoted as ((.)). For instance, /2 is estimated as 

/2=.((-^£cos(2(^,))), (2.3) 

where k labels the particles the event. The equality becomes strict as the number of events ap- 
proaches infinity, which is assumed implicitly. In the variable-axes case we rotate the particles by 
the angle 0* in each event. Thus 

/2* = £* = ((-i:cos[2((/>,-r)])). (2.4) 

The rotation angle 0* depends itself on the distribution of particles in the given event. By definition, 
it is chosen in such a way that the quantity ^ L^^j cos[2(^<: — 0*)] assumes maximum, which gives 
the conditions 



co^{l^*) = Y2/^Yi+Xl ^m{l^*)=X2/\lY^+Xl (2.5) 

Y ^ 1 " 

72 = - V cos(2(/)<.), ^2 = - V sin(20i). 



Using the above formulas in Eq. ( p.4| ) yields 



n = (('Jn+^i)) = (( 



\ 2 / \ 2 

1 " , \ /i « 



£ cos(20,) + - £ sin(20,) )). (2.6) 



We see that the variable-axes moment corresponds to an average of the square root of sums ([2.5|), 
thus is a highly "non-local" object, involving upon expansion infinitely many fixed-axes moments. 

For sufficiently large n one may evaluate Eq. ( ^^ with the help of the central limit theorem. 
Consider the variables Ck = cos{l^t) and Sk = sin{l(j)k). Their averages and variances are 

c = e, o^c=\- e^ ^ = 0, = ^. (2.7) 

Importantly, there is no correlation between Y2 and X2, as ^ /q^'^ <i0 cos(2(^) sin(20)/(0) = 0. 
Thus, according to the central limit theorem, the distribution of Y2 and X2 is Gaussian. Introducing 

72=^cosa, X2 = ^sina, q^ = Y^+Xl, g = ^ - ^ = ^ ^^^^ - 1, (2.8) 
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^2 n = 10, 20, 50, 100, 400 CJ ( f |) / f 2 n = 10, 20, 50, 100, 400 




IkI^ otl ots ^ o!l 0^2 0^3 0^4 ots ^ 



Figure 2: Toy model. Left: Dependence of the variable-axes moment /| on the fixed-axes quadrupole 
moment e for several values of the number of particles n. As n increases, we pass from top to bottom with 
the presented curves. The straight line is the « — > oo limit, i.e. /| = e. Right: the same for the scaled standard 
deviation. 



we may write this distribution in the form 



f{X2,Y2)=fiq,a) 



7rVl-2£2 



exp 



+ E'^ — Iqecosa 



+ n5q^ sin^ a 



(2.9) 



We need below the integral of this distribution over a, which can be expanded as \ \13[ |11| ] 



27t 



daf{q,a)- 



2n 



7rVl-2£2 



exp 



q^ + £^ 
l-2£2 



2n£q 



, (2.10) 



where Ij{x) are the modified Bessel functions. We may now express Eq. (2.4) as the series involving 
the confluent hypergeometric function, 



/I = qdqdaqf{q,a)-- 



l-2£''^,.„2^i^(J + i)^(7 + |) 



I (2s 



7! 



|2 



, (2.11) 



which converges fast and can be used for practical calculations in a truncated form. At £ = 
(azimuthally symmetric distribution) we have the very simple result 



(2.12) 



which shows the expected 1 / y/n behavior for a statistical fluctuation. The numerical results ob- 
tained with the series ( 2.11 ) are presented in Fig. left side. We note that the effect of the departure 
of /I from £ is strongest at low £ and low n. 

The evaluation of the second moment in the q variable yields 

1 + (?1-1)£2 



J qdqdaq^f{q,a) 



(2.13) 



From Eqs. (2.11,2.13) we can now obtain the variance of the distribution of the variable-axes mo- 
ment. Again, a simple formula follows for the case £ = 0, where var(/|) = (1 — f )/«• The scaled 
variance and scaled standard deviation are 

var(/2*) _ 7^-— oif*) _ 



/I 



/I 



n 



1 -0.523, (£ = 0). 



(2.14) 
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Note that in this case there is no dependence of the scaled standard deviation on n. The case of 
general e obtained numerically for various values of n is shown in Fig. right side. According to 



Eq. ( |2.14| ), all curves approach the limit — 1 as £ ^ 0. At the other end, in the limit of ne^ 



we have the expansions /| = £ + 1/ {4en) + ... and a(/2*)//2* = [1/ (2£) -e]/n + .... 
3. The general case 

In the general case the analysis can be carried out in full analogy to the toy model For 
simplicity, in our analytic study we neglect correlations between locations of sources. If such corre- 
lations are strong, their analytic inclusion is difficult and one has to resort to numerical simulations 
such as those presented below. Compared to the toy model, the full two-dimensional case involves 
the fixed-axes moments 4./ = Jo'^^pdpfi{p)p'^ /n, where n is the number of sources. We have 



chosen the normalization J^2Kpdpfo{p) = n. Generally, in analogy to Eq. ( 2.11 ) 



V2a, 



2 ~ r(m+i)r(m + |)iFi(-i;m+l;-^) 



.-^ ^ I (2502)- iA, (3.1) 

where 

2^2 2 1 11 

^2=4,2, o-y^ = ^feo -24,2 +^2m)' '^X2 = i^ihkfi - hkA) , ^^2o^ ~2o^' 

For the special case of central collisions we have the very simple results 



s* = #P, ^ = ./^c. 0.523, (^ = 0). (3.3) 

Since correlations between the location of sources effectively reduce the number of sources n, they 
lead to an increase of £*, but keep its ^-independent scaled variance practically constant, as shown 
by the simulations of the next Section. Ref. [|l|] contains more discussion. 

4. Numerical simulations in various Glauber models 

We have studied a few variants of Glauber-like models. In the standard wounded nucleon 
model [^] the weight w = 1/2 is attributed to the point in the transverse plane at the position of the 
wounded nucleon. The wounding cross section is 42 mb. For binary collisions the weight w = 1 is 
attributed to each collision point. We remark that only relative magnitude of weights is important in 
studies of fluctuations. A successful description of the RHIC multiplicities has been achieved with 
a mixed model, amending wounded nucleons with some admixture of binary collisions p^ ]. 
Then the wounded nucleon obtains the weight w = (1 — a)/2, and the binary collision the weight 
w = a. The total weight averaged over events is then (1 — a)N^/2 + aA'bin- The fits to particle 
multiplicities of Ref. |[l^] give a = 0. 145 at Js^n = 200 GeV. We also consider a model with 



hot spots in the spirit of Ref. [18], assuming that the cross section for a semi-hard binary colli- 
sions producing a hot-spot is tiny, ahot-spot = 0.5 mb, however when such a rare collision occurs 
it produces on the average a very large amount of the transverse energy equal to aaw/ Chot-spot- 
Each source from the previously described models deposits the transverse energy with a certain 
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Figure 3: Left: The moment e and its scaled standard deviation for the analyzed models plotted as functions 
of the impact parameter. Gold-gold collisions. Right: the same for e*. The results for the color-glass 
condensate come from Ref. | ]I4| ] (dot-dashed line in top left figure) and Ref. [ pO| ] (circles in lower right 
figure). 



probability distribution. To incorporate this effect, we superimpose the F distribution, multiplying 
the weights of the considered model with the randomly distributed number from the gamma distri- 
bution g{w, k) = w"^^^ K"'^exp(— K"w)/r(K"). Here we do this superposition on the hot-spot model, 
labeled hot-spot+T. Thus, we take the weights (1 — a)g{w, K")/2 for the wounded nucleons and 
(Xg{w, K")aw/cJhot-spot for the binary colhsions. We set k = 0.5, which gives var(w) = 5. The four 
considered models (wounded nucleon, mixed, hot-spot, and hot-spot-i-F) differ substantially by the 
number of sources and the amount of the built-in fluctuations. 

We observe that in all four models e is practically independent of the model (top left panel 
of Fig. On the other hand, the scaled standard deviation (lower left panel of Fig. ||) displays 
a strong dependence on the model at low values of b, with the hot-spot-i-F model yielding about 
twice as much as the mixed model. We also notice a very strong dependence onb. At = the 
curves diverge due to dividing by the vanishing value of £. The fluctuations are larger in models 
effectively having the lower number of sources, which is obvious from the statistical point of view. 

As already noted in Refs. |14, 19], the value of e obtained with the color glass condensate 
(CGC) is substantially higher than in all Glauber-like models reported here (upper curve in the left 
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top panel of Fig. |3|). 

The quadrupole moment e* and its scaled standard deviation are show on the right side of Fig. || 
We observe a strong model dependence of e* at low values of b, with models having effectively 
lower number of sources yielding higher values. Atb = the hot-spot+F model yields three times 
more than the wounded nucleon model. For all models the scaled standard deviation is close to 



the value 0.5 for central collisions (in agreement with the results (|3.3|)) and drops to about 0.3 
at b = 14 fm. At intermediate values of b the relative difference in a(£*)/e* between various 
considered models is at the level of 10-15%, which is not a very strong effect. The CGC result of 
Ref. [20] is lower than in the Glauber models (circles in the right bottom panel of Fig. 
The harmonic profiles //(p) and //(p) are displayed in Ref. [|I|] 



5. Jet quenching 

We have used the model of Refs. [ ^T| , ^ of the jet energy loss in order to explore the role 
of the event-by-event rotated absorbing medium. In order to take into account the variable-axes 
geometry, we use /*(p,0) as the density of the scattering centers for the propagating parton. The 
rare jet production event is distributed according to the fixed-axes profile, see Fig. left side. The 
resulting increase of the eccentricity of the absorbent is expected to increase the asymmetry of the 
jet absorption. A very similar effect has been discussed for profiles calculated in the CGC model 
[p3[], where an increase in V2 by about 10 — 15% has been found. The absorbing medium formed in 
each event is rotated and also shifted. The elliptic flow (see Fig. |^ right side) at centralities larger 
than 20% calculated with the wounded-nucleon model in the fixed-axes frame (solid line), which 
serves as a reference, comes out similar to the result of the hot-spot model in the variable-axes 
frame (dashed line). Only if the shift and rotation of the opaque medium were neglected (dotted 
line) the modification of the shape leads to an increase of the high pj elliptic flow coefficient V2 by 
about 10— 15%. The cancellation of the effects of the increased eccentricity of the medium and 



absorbent 




,0.1 - 



0.05 - 






400 



Figure 4: Left: the variable axes geometry of the absorbent vs. the fixed-axes geometry of the jet production 
point. Right: V2 at high pr as a function of the number of wounded nucleons, obtained with the variable- 
axes density /* (p ) for the hot-spot scenario (dashed line), and with the fixed-axes density of the wounded 
nucleons /(p) (solid line). The dotted line represents the result for the variable-axes density but without the 
shift and rotation of the opaque medium. 



7 



Fluctuations from Glauber models 



Wojciech Broniowski 



of the shift and rotation happens also for the other considered models (at larger centralities). The 
rotation of the absorbing medium yields about 2/3, and the shift about 1/3 of the total cancellation 
effect. 



6. Fluctuations of the elliptic flow 



The fluctuations of the elliptic flow, which are an important probe of the nature of the early- 
stage dynamics of the system [24], have recently been measured at RHIC [ |I0| , [Tl], 12]. The ex- 
perimental procedure used in these analyses identifies the elliptic flow coefficient with the variable 
axes V2, here denoted as Vj. The relevance of studies of fluctuations of the initial shape comes from 
the well-known fact that for small elliptic asymmetry one expects on hydrodynamic grounds the 
relation 

^ = il(p. (6.1, 



As argued in Ref. [|25|], the result ( p.lD indicates that the mean free path in the matter created in 
the initial stages of the heavy-ion collisions is very small, although turbulence does not develop. 
Comparison of the data to our Glauber calculations is made in Fig. |[ For central collisions we 
expect 



:0) ~ 



1-0.52, 



(6.2) 



which is compatible to the data, although the error bars are large. 

We end this talk with some comments on the derivation of Eq. (5J) as well as on higher- 
harmonic probes. Perturbation theory applied to smooth {i.e. linearization is sensible) hydrody- 
namics together with hierarchy of relaxation times for subsequent harmonics, e.g. I2 3> T4, leads 
to further results [|l|]. In particular. 



,,* „*2 ,,*2 
V4 ~ £ ~ V2 . 



(6.3) 



In Ref. the variable V4 /v\ has been suggested as a sensitive probe of the hydrodynamic evo- 



lution. The simulations of Refs. [ ]26| , [27[ ] show that with increasing time the value of V2 saturates. 




■PHOBOS 
STAR 
wounded AsVs* 

hot spot + r A8*/s* 



50 100 150 200 250 300 350 m 

Figure 5: Fluctuations of v^. Data from Refs. |ri[ |l2[] . 
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while V4 quickly assumes the value proportional to V2, supporting the assumption T2 » T4 used in 



the above argumentation. For the fluctuations one gets immediately from Eq. ( |6.3| ) the prediction 



^.2^. (6.4, 



Relation (5.4), if verified experimentally, would support the scenario of smooth hydro evolution 
with the mentioned hierarchy of scales. On similar grounds, for the azimuthal Hanbury-Brown- 
Twiss (HBT) correlation radius, /?™^((/)), one expects 

where = R™'^ + 2/?^"^^ cos (20) + 2/?^*'^'^ cos (40) + . . . . 

7. Conclusion 

Here are our main points: 

• We have analyzed four Glauber-like models, with different degree of fluctuation: the wounded- 
nucleon model, the mixed model, the hot-spot model, and the hot-spot model with the super- 
imposed r distribution. 

• We have obtained numerically the fixed-axes and variable-axes harmonic profiles ^ and 
analyzed their moments. The variable-axes moments £*, and the fixed-axes scaled standard 
deviation o{e)/e are sensitive to the choice of the model, while o{e*)/e* is not, changing 
at most by 10-15% from model to model at intermediate values of b. 

• Analytic formulas explain certain features of the simulations, in particular, they show that at 
b = the multiple-axes scaled variances are close to the value 0.5, insensitive of the model 
used, the collision energy, the mass number of the colliding nuclei, or the number of particle 
sources. The behavior of a (£*)/£* at low b is thus largely governed by the statistics 

• Fixed-reaction-plane experimental analyses would reveal more information on the system 
and would allow to discriminate the theoretical predictions, as fluctuations of e are sensitive 
to the chosen model. 

• For the jet V2 we find that the effect of the increased variable-axes eccentricity is largely 
canceled by the shift of the center of mass and the rotation of the principal axes of the 
absorbing medium. This leads to practically no change of the jet emission asymmetry at 
intermediate and large impact parameters. Only at small b the increase of the quadrupole 
moment takes over the relatively less important shift and rotation. 

• On hydrodynamic grounds, the analysis of the variable-axes moments in the coordinate space 



carries over to the collective flow and analysis of v^- In particular, Eq. ( |6.2[ ) holds for the 
variable-axes elliptic flow coefficient. 

Under assumptions of smoothness, perturbation theory made on top of azimuthally symmet- 
ric hydro leads to sensitivity of higher-harmonic late-time measures, v|, etc., to the 
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initial quadrupole deformation £*{to) only. Higher harmonics of the initial shape deforma- 
tion are irrelevant, as they presumably are damped fast. A number of relations follows for 
various measures and their event-by-event fluctuations, e.g. Eq. ( |6.4[ ). 



It would be a challenge to measure the v\ fluctuations and test the smooth hydro assumption 



by verifying relation ( |6.4[ ). 



One of us (WB) thanks Paul Sorensen, Constantin Loizides, and Wit Busza for useful discus- 
sions concerning the experimental determination of V2 and its fluctuations. 
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